Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria).

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

Phagocytosis underpins the biotrophic lifestyle of intracellular parasites in the class Phytomyxea (Rhizaria).

Phytomyxea are intracellular biotrophic parasites infecting plants and stramenopiles, including the agriculturally impactful Plasmodiophora brassicae and the brown seaweed pathogen Maullinia ectocarpii. They belong to the clade Rhizaria, where phagotrophy is the main mode of nutrition. Phagocytosis is a complex trait of eukaryotes, well documented for free-living unicellular eukaryotes and specific cellular types of animals. Data on phagocytosis in intracellular, biotrophic parasites are scant. Phagocytosis, where parts of the host cell are consumed at once, is seemingly at odds with intracellular biotrophy. Here we provide evidence that phagotrophy is part of the nutritional strategy of Phytomyxea, using morphological and genetic data (including a novel transcriptome of M. ectocarpii). We document intracellular phagocytosis in P. brassicae and M. ectocarpii by transmission electron microscopy and fluorescent in situ hybridization. Our investigations confirm molecular signatures of phagocytosis in Phytomyxea and hint at a small specialized subset of genes used for intracellular phagocytosis. Microscopic evidence confirms the existence of intracellular phagocytosis, which in Phytomyxea targets primarily host organelles. Phagocytosis seems to coexist with the manipulation of host physiology typical of biotrophic interactions. Our findings resolve long debated questions on the feeding behaviour of Phytomyxea, suggesting an unrecognized role for phagocytosis in biotrophic interactions.

A Mitosome With Distinct Metabolism in the Uncultured Protist Parasite Paramikrocytos canceri (Rhizaria, Ascetosporea).

Ascetosporea are endoparasites of marine invertebrates that include economically important pathogens of aquaculture species. Owing to their often-minuscule cell sizes, strict intracellular lifestyle, lack of cultured representatives and minimal availability of molecular data, these unicellular parasites remain poorly studied. Here, we sequenced and assembled the genome and transcriptome of Paramikrocytos canceri, an endoparasite isolated from the European edible crab Cancer pagurus. Using bioinformatic predictions, we show that P. canceri likely possesses a mitochondrion-related organelle (MRO) with highly reduced metabolism, resembling the mitosomes of other parasites but with key differences. Like other mitosomes, this MRO is predicted to have reduced metabolic capacity and lack an organellar genome and function in iron-sulfur cluster (ISC) pathway-mediated Fe-S cluster biosynthesis. However, the MRO in P. canceri is uniquely predicted to produce ATP via a partial glycolytic pathway and synthesize phospholipids de novo through the CDP-DAG pathway. Heterologous gene expression confirmed that proteins from the ISC and CDP-DAG pathways retain mitochondrial targeting sequences that are recognized by yeast mitochondria. This represents a unique combination of metabolic pathways in an MRO, including the first reported case of a mitosome-like organelle able to synthesize phospholipids de novo. Some of these phospholipids, such as phosphatidylserine, are vital in other protist endoparasites that invade their host through apoptotic mimicry.

Morphological and molecular diversity of monothalamids (Rhizaria, Foraminifera), including two new species and a new genus, from SW Greenland.

Single-chambered (monothalamous) foraminifera are poorly known compared to their multichambered relatives. In this first study of monothalamids from Greenland, we describe one new genus and two new species belonging to different clades from the Nuuk fjord system. Nujappikia idaliae Gooday & Holzmann gen. nov. sp. nov. (Clade Y) has a bottle-shaped test terminating in a single aperture located on a short neck. The flexible wall is basically organic but with a very fine agglutinated veneer. Bathyallogromia kalaallita Gooday & Holzmann sp. nov. (Clade C) has a broadly ovate test with an organic wall and a mound-like apertural structure. It is larger and genetically distinct from the two other Bathyallogromia species, both from the Southern Ocean. A survey of the morphological diversity of monothalamids in our samples revealed 49 morphospecies, of which 19, including the two new species, yielded DNA sequences. Five were assigned to the genera Bathysiphon, (Clade BM), Micrometula. (Clade BM), Psammophaga. (Clade E), Hippocrepinella (Clade D) and Crithionina (Clade J). The remaining twelve represented unknown taxa branching in clades A, C, F, and Y and one new clade. Our results add to growing evidence that monothalamids are common and diverse in fjords and other high-latitude settings.

Taxon-rich transcriptomics supports higher-level phylogeny and major evolutionary trends in Foraminifera.

Foraminifera, classified in the supergroup Rhizaria, are a common and highly diverse group of mainly marine protists. Despite their evolutionary and ecological importance, only limited genomic data (one partial genome and nine transcriptomic datasets) have been published for this group. Foraminiferal molecular phylogeny is largely based on 18S rRNA gene sequence analysis. However, due to highly variable evolutionary rates of substitution in ribosomal genes plus the existence of intragenomic variation at this locus, the relationships between and within foraminiferal classes remain uncertain. We analyze transcriptomic data from 28 species, adding 19 new species to the previously published dataset, including members of the strongly under-represented class Monothalamea. A phylogenomic reconstruction of Rhizaria, rooted with alveolates and stramenopiles, based on 199 genes and 68 species supports the monophyly of Foraminifera and their sister relationship to Polycystinea. The phylogenomic tree of Foraminifera is very similar to the 18S rRNA tree, with the paraphyletic single-chambered monothalamids giving rise to the multi-chambered Tubothalamea and Globothalamea. Within the Monothalamea, our analyses confirm the monophyly of the giant, deep-sea xenophyophores that branch within clade C and indicate the basal position of monothalamous clades D and E. The multi-chambered Globothalamea are monophyletic and comprise the paraphyletic Textulariida and monophyletic Rotaliida. Our phylogenomic analyses support major evolutionary trends of Foraminifera revealed by ribosomal phylogenies and reinforce their current higher-level classification.

Intra-genomic rRNA gene variability of Nassellaria and Spumellaria (Rhizaria, Radiolaria) assessed by Sanger, MinION and Illumina sequencing.

Ribosomal RNA (rRNA) genes are known to be valuable markers for the barcoding of eukaryotic life and its phylogenetic classification at various taxonomic levels. The large-scale exploration of environmental microbial diversity through metabarcoding approaches has been focused mainly on the V4 and V9 regions of the 18S rRNA gene. The accurate interpretation of such environmental surveys is hampered by technical (e.g. PCR and sequencing errors) and biological biases (e.g. intra-genomic variability). Here we explored the intra-genomic diversity of Nassellaria and Spumellaria specimens (Radiolaria) by comparing Sanger sequencing with Illumina and Oxford Nanopore Technologies (MinION). Our analysis determined that intra-genomic variability of Nassellaria and Spumellaria is generally low, yet some Spumellaria specimens showed two different copies of the V4 with <97% similarity. Of the different sequencing methods, Illumina showed the highest number of contaminations (i.e. environmental DNA, cross-contamination, tag-jumping), revealed by its high sequencing depth; and MinION showed the highest sequencing rate error (~14%). Yet the long reads produced by MinION (~2900 bp) allowed accurate phylogenetic reconstruction studies. These results highlight the requirement for a careful interpretation of Illumina-based metabarcoding studies, in particular regarding low abundant amplicons, and open future perspectives towards full-length rDNA environmental metabarcoding surveys.

Phylogenetic analysis confirms the taxonomic placement of the marine flagellate Hermesinum adriaticum (Thecofilosea, Cercozoa, Rhizaria).

Hermesinum adriaticum is a rare marine and brackish flagellate that is of considerable interest due to its markable and fossilizable siliceous skeleton. Based on this skeleton, Hermesinum was initially considered a microalga of the Dictyochophyceae (Ochrophyta, Stramenopiles). Later on, it was assigned to the Ebriida due to its similarity to Ebria tripartita. The taxonomic assignment of the Ebriida, however, changed several times until it was placed within the Thecofilosea (Cercozoa, Rhizaria), based on genetic data of E. tripartita. We sequenced the 18S marker gene sequence of Hermesinum and confirm the close relationship of Ebria and Hermesinum.

Transfer of the thecate amoebae Lecythium spinosum and Pamphagus armatus to Rhizaspis (Thecofilosea, Cercozoa, Rhizaria).

Thecofilosea is a class in Cercozoa (Rhizaria) comprising mainly freshwater-inhabiting algivores. Recently, numerous isolates of thecofilosean amoebae have been cultured and were characterized by an integrated morphological and molecular approach. The captivating spine-bearing taxa in Thecofilosea were not yet molecularly characterized due to being very rare. There are only two known spine-bearing species, Pamphagus armatus and Lecythium spinosum, which were synonymized by Penard in 1902. Due to a morphological difference of those taxa, we discuss here that we disagree with this taxonomical act. We further isolated single cells of Pamphagus armatus directly from their habitat and successfully sequenced their SSU rDNA, which we subjected to phylogenetic analyses. We show that Pamphagus armatus branches within the Rhizaspididae (Tectofilosida, Thecofilosea). Accordingly, we transfer Pamphagus armatus and the assumingly closely related species Lecythium spinosum to Rhizaspis.

A Morpho-molecular Perspective on the Diversity and Evolution of Spumellaria (Radiolaria).

Spumellaria (Radiolaria, Rhizaria) are holoplanktonic amoeboid protists, ubiquitous and abundant in the global ocean. Their silicified skeleton preserves very well in sediments, displaying an excellent fossil record extremely valuable for paleo-environmental reconstruction studies, from where most of their extant diversity and ecology have been inferred. This study represents a comprehensive classification of Spumellaria based on the combination of ribosomal taxonomic marker genes (rDNA) and morphological characteristics. In contrast to established taxonomic knowledge, we demonstrate that symmetry of the skeleton takes more importance than internal structures at high classification ranks. Such reconsideration allows gathering different morphologies with concentric structure and spherical or radial symmetry believed to belong to other Radiolaria orders from the fossil record, as for some Entactinaria families. Our calibrated molecular clock dates the origin of Spumellaria in the middle Cambrian (ca. 515 Ma), among the first radiolarian representatives in the fossil record. This study allows a direct connection between living specimens and extinct morphologies from the Cambrian, bringing both a standpoint for future molecular environmental surveys and a better understanding for paleo-environmental reconstruction analysis.

New Permian radiolarians from east Asia and the quantitative reconstruction of their evolutionary and ecological significances.

The biostratigraphically important Permian radiolarian genera Pseudoalbaillella sensu stricto and Follicucullus (Follicucullidae, Polycystinea) are discriminated by morphological gaps in their wings and segmentation. Previous statistical analyses demonstrated that Longtanella fills morphological gaps between these two genera. Longtanella has long been regarded as a junior synonym of Parafollicucullus, and only a few species have been described. Herein several true Longtanella species are recognized from South China, and eight new species and five indeterminate species are described and illustrated to prove the validity of the genus Longtanella. In addition, a new genus, Parafollicucullinoides gen. nov., is described. Their palaeogeographic distributions and living environments are explored by applying correspondence analysis (CA), with occurrence datasets of selected fusulinacean genera from the Japanese Islands, China and Sundaland. CA results indicate that Longtanella was present to a limited extent in warmer conditions in the fusulinacean Province B and C during Kungurian-Roadian time, and possibly lived above the thermocline and below the deepest limit of fusulinaceans. The Pseudoalbaillella and the Follicucullus group preferred open ocean conditions, living below the thermocline and distributed not only in the 'Equatorial Warm Water Province', but also the northern peri-Gondwana Cool Water Province and the southern North Cool Water Province.

Microbial dynamics of elevated carbon flux in the open ocean's abyss.

In the open ocean, elevated carbon flux (ECF) events increase the delivery of particulate carbon from surface waters to the seafloor by severalfold compared to other times of year. Since microbes play central roles in primary production and sinking particle formation, they contribute greatly to carbon export to the deep sea. Few studies, however, have quantitatively linked ECF events with the specific microbial assemblages that drive them. Here, we identify key microbial taxa and functional traits on deep-sea sinking particles that correlate positively with ECF events. Microbes enriched on sinking particles in summer ECF events included symbiotic and free-living diazotrophic cyanobacteria, rhizosolenid diatoms, phototrophic and heterotrophic protists, and photoheterotrophic and copiotrophic bacteria. Particle-attached bacteria reaching the abyss during summer ECF events encoded metabolic pathways reflecting their surface water origins, including oxygenic and aerobic anoxygenic photosynthesis, nitrogen fixation, and proteorhodopsin-based photoheterotrophy. The abundances of some deep-sea bacteria also correlated positively with summer ECF events, suggesting rapid bathypelagic responses to elevated organic matter inputs. Biota enriched on sinking particles during a spring ECF event were distinct from those found in summer, and included rhizaria, copepods, fungi, and different bacterial taxa. At other times over our 3-y study, mid- and deep-water particle colonization, predation, degradation, and repackaging (by deep-sea bacteria, protists, and animals) appeared to shape the biotic composition of particles reaching the abyss. Our analyses reveal key microbial players and biological processes involved in particle formation, rapid export, and consumption, that may influence the ocean's biological pump and help sustain deep-sea ecosystems.

Shifts in the protist community associated with an anticyclonic gyre in the Alboran Sea (Mediterranean Sea).

The diversity of protists was researched in the Alboran Sea (SW Mediterranean Sea) by means of high-throughput sequencing technologies based on the amplification of the V9 region of 18S rRNA. Samples were collected at different depths in seven stations following an environmental gradient from a coastal upwelling zone to the core of an oligotrophic anticyclonic gyre (AG). Sampling was performed during summer, when the water column was stratified. The superphyla Alveolata, Stramenopila and Rhizaria accounted for 84% of the total operational taxonomic units (OTUs). The most diverse groups were Dinophyceae (21% of OTUs), Marine Alveolates-II (MALV-II; 20%), Ciliophora (9%) and MALV-I (6%). In terms of read abundance, the predominant groups were Dinophyceae (29%), Bacillariophyta (14%), MALV-II (11%) and Ciliophora (11%). Samples were clustered into three groups according to the sampling depth and position. The shallow community in coastal stations presented distinguishable patterns of diatoms and ciliates compared with AG stations. These results indicate that there was a strong horizontal coupling between phytoplankton and ciliate communities. Abundance of Radiolaria and Syndiniales increased with depth. Our analyses demonstrate that the stratification disruption produced by the AG caused shifts in the trophic ecology of the plankton assemblages inducing a transition from bottom-up to top-down control.

Assessing the response of micro-eukaryotic diversity to the Great Acceleration using lake sedimentary DNA.

Long-term time series have provided evidence that anthropogenic pressures can threaten lakes. Yet it remains unclear how and the extent to which lake biodiversity has changed during the Anthropocene, in particular for microbes. Here, we used DNA preserved in sediments to compare modern micro-eukaryotic communities with those from the end of the 19th century, i.e., before acceleration of the human imprint on ecosystems. Our results obtained for 48 lakes indicate drastic changes in the composition of microbial communities, coupled with a homogenization of their diversity between lakes. Remote high elevation lakes were globally less impacted than lowland lakes affected by local human activity. All functional groups (micro-algae, parasites, saprotrophs and consumers) underwent significant changes in diversity. However, we show that the effects of anthropogenic changes have benefited in particular phototrophic and mixotrophic species, which is consistent with the hypothesis of a global increase of primary productivity in lakes.

Long transposon-rich centromeres in an oomycete reveal divergence of centromere features in Stramenopila-Alveolata-Rhizaria lineages.

Centromeres are chromosomal regions that serve as platforms for kinetochore assembly and spindle attachments, ensuring accurate chromosome segregation during cell division. Despite functional conservation, centromere DNA sequences are diverse and often repetitive, making them challenging to assemble and identify. Here, we describe centromeres in an oomycete Phytophthora sojae by combining long-read sequencing-based genome assembly and chromatin immunoprecipitation for the centromeric histone CENP-A followed by high-throughput sequencing (ChIP-seq). P. sojae centromeres cluster at a single focus at different life stages and during nuclear division. We report an improved genome assembly of the P. sojae reference strain, which enabled identification of 15 enriched CENP-A binding regions as putative centromeres. By focusing on a subset of these regions, we demonstrate that centromeres in P. sojae are regional, spanning 211 to 356 kb. Most of these regions are transposon-rich, poorly transcribed, and lack the histone modification H3K4me2 but are embedded within regions with the heterochromatin marks H3K9me3 and H3K27me3. Strikingly, we discovered a Copia-like transposon (CoLT) that is highly enriched in the CENP-A chromatin. Similar clustered elements are also found in oomycete relatives of P. sojae, and may be applied as a criterion for prediction of oomycete centromeres. This work reveals a divergence of centromere features in oomycetes as compared to other organisms in the Stramenopila-Alveolata-Rhizaria (SAR) supergroup including diatoms and Plasmodium falciparum that have relatively short and simple regional centromeres. Identification of P. sojae centromeres in turn also advances the genome assembly.

Making sense of environmental sequencing data: Ecologically important functional traits of the protistan groups Cercozoa and Endomyxa (Rhizaria).

We have compiled a database of functional traits for two widespread and ecologically important groups of protists, Cercozoa and Endomyxa (Rhizaria). The functional traits of microorganisms are crucially important for interpreting results from environmental sequencing surveys. Linking morphological and ecological traits to environmental factors is common practice in studies involving micro- and macroorganisms, but is rarely applied to protists. Our database provides functional and ecologically significant traits linked to morphology, nutrition, locomotion and habitats. We discuss how the use of functional traits may help to unveil underlying ecosystem processes. This database is intended as a common reference for the molecular ecology community and will boost the understanding of ecosystem functions, especially those driven by biological interactions.

Time Calibrated Morpho-molecular Classification of Nassellaria (Radiolaria).

Nassellaria are marine protists belonging to the Radiolaria lineage (Rhizaria). Their skeleton, made of opaline silica, exhibit an excellent fossil record, extremely valuable in micro-paleontological studies for paleo-environmental reconstruction. Yet, to date very little is known about the extant diversity and ecology of Nassellaria in contemporary oceans, and most of it is inferred from their fossil record. Here we present an integrative classification of Nassellaria based on taxonomical marker genes (18S and 28S ribosomal DNA) and morphological characteristics obtained by optical and scanning electron microscopy imaging. Our phylogenetic analyses distinguished 11 main morpho-molecular clades relying essentially on the overall morphology of the skeleton and not on internal structures as previously considered. Using fossil calibrated molecular clock we estimated the origin of Nassellaria among radiolarians primitive forms in the Devonian (ca. 420 Ma), that gave rise to living nassellarian groups in the Triassic (ca. 250 Ma), during the biggest diversification event over their evolutionary history. This morpho-molecular framework provides both a new morphological classification easier to identify under light microscopy and the basis for future molecular ecology surveys. Altogether, it brings a new standpoint to improve our scarce understanding of the ecology and worldwide distribution of extant nassellarians.

High contribution of Rhizaria (Radiolaria) to vertical export in the California Current Ecosystem revealed by DNA metabarcoding.

Passive sinking of particulate organic matter (POM) is the main mechanism through which the biological pump transports surface primary production to the ocean interior. However, the contribution and variability of different biological sources to vertical export is not fully understood. Here, we use DNA metabarcoding of the 18S rRNA gene and particle interceptor traps (PITs) to characterize the taxonomic composition of particles sinking out of the photic layer in the California Current Ecosystem (CCE), a productive system with high export potential. The PITs included formalin-fixed and 'live' traps to investigate eukaryotic communities involved in the export and remineralization of sinking particles. Sequences affiliated with Radiolaria dominated the eukaryotic assemblage in fixed traps (90%), with Dinophyta and Metazoa making minor contributions. The prominence of Radiolaria decreased drastically in live traps, possibly due to selective consumption by copepods, heterotrophic nanoflagellates, and phaeodarians that were heavily enriched in these traps. These patterns were consistent across the water masses surveyed extending from the coast to offshore, despite major differences in productivity and trophic structure of the epipelagic plankton community. Our findings identify Radiolaria as major actors in export fluxes in the CCE.

Spatial Variability of Picoeukaryotic Communities in the Mariana Trench.

Picoeukaryotes play prominent roles in the biogeochemical cycles in marine ecosystems. However, their molecular diversity studies have been confined in marine surface waters or shallow coastal sediments. Here, we investigated the diversity and metabolic activity of picoeukaryotic communities at depths ranging from the surface to the abyssopelagic zone in the western Pacific Ocean above the north and south slopes of the Mariana Trench. This was achieved by amplifying and sequencing the V4 region of both 18S ribosomal DNA and cDNA using Illumina HiSeq sequencing. Our study revealed: (1) Four super-groups (i.e., Alveolata, Opisthokonta, Rhizaria and Stramenopiles) dominated the picoeukaryote assemblages through the water column, although they accounted for different proportions at DNA and cDNA levels. Our data expand the deep-sea assemblages from current bathypelagic to abyssopelagic zones. (2) Using the cDNA-DNA ratio as a proxy of relative metabolic activity, the highest activity for most subgroups was usually found in the mesopelagic zone; and (3) Population shift along the vertical scale was more prominent than that on the horizontal differences, which might be explained by the sharp physicochemical gradients along the water depths. Overall, our study provides a better understanding of the diversity and metabolic activity of picoeukaryotes in water columns of the deep ocean in response to varying environmental conditions.

Multigene phylogeny and cell evolution of chromist infrakingdom Rhizaria: contrasting cell organisation of sister phyla Cercozoa and Retaria.

Infrakingdom Rhizaria is one of four major subgroups with distinct cell body plans that comprise eukaryotic kingdom Chromista. Unlike other chromists, Rhizaria are mostly heterotrophic flagellates, amoebae or amoeboflagellates, commonly with reticulose (net-like) or filose (thread-like) feeding pseudopodia; uniquely for eukaryotes, cilia have proximal ciliary transition-zone hub-lattices. They comprise predominantly flagellate phylum Cercozoa and reticulopodial phylum Retaria, whose exact phylogenetic relationship has been uncertain. Given even less clear relationships amongst cercozoan classes, we sequenced partial transcriptomes of seven Cercozoa representing five classes and endomyxan retarian Filoreta marina to establish 187-gene multiprotein phylogenies. Ectoreta (retarian infraphyla Foraminifera, Radiozoa) branch within classical Cercozoa as sister to reticulose Endomyxa. This supports recent transfer of subphylum Endomyxa from Cercozoa to Retaria alongside subphylum Ectoreta which embraces classical retarians where capsules or tests subdivide cells into organelle-containing endoplasm and anastomosing pseudopodial net-like ectoplasm. Cercozoa are more homogeneously filose, often with filose pseudopodia and/or posterior ciliary gliding motility: zooflagellate Helkesimastix and amoeboid Guttulinopsis form a strongly supported clade, order Helkesida. Cercomonads are polyphyletic (Cercomonadida sister to glissomonads; Paracercomonadida deeper). Thecofilosea are a clade, whereas Imbricatea may not be; Sarcomonadea may be paraphyletic. Helkesea and Metromonadea are successively deeper outgroups within cercozoan subphylum Monadofilosa; subphylum Reticulofilosa (paraphyletic on site-heterogeneous trees) branches earliest, Granofilosea before Chlorarachnea. Our multiprotein trees confirm that Rhizaria are sisters of infrakingdom Halvaria (Alveolata, Heterokonta) within chromist subkingdom Harosa (= SAR); they further support holophyly of chromist subkingdom Hacrobia, and are consistent with holophyly of Chromista as sister of kingdom Plantae. Site-heterogeneous rDNA trees group Kraken with environmental DNA clade 'eSarcomonad', not Paracercomonadida. Ectoretan fossil dates evidence ultrarapid episodic stem sequence evolution. We discuss early rhizarian cell evolution and multigene tree coevolutionary patterns, gene-paralogue evidence for chromist monophyly, and integrate this with fossil evidence for the age of Rhizaria and eukaryote cells, and revise rhizarian classification.

Rhizarian 'Novel Clade 10' Revealed as Abundant and Diverse Planktonic and Terrestrial Flagellates, including Aquavolon n. gen.

Rhizarian 'Novel Clade 10' (NC10) is frequently detected by 18S rRNA gene sequencing studies in freshwater planktonic samples. We describe a new genus and two species of eukaryovorous biflagellate protists, Aquavolon hoantrani n. gen. n. sp. and A. dientrani n. gen. n. sp., which represent the first morphologically characterized members of NC10, here named Aquavolonida ord. nov. The slightly metabolic cells possess naked heterodynamic flagella, whose kinetosomes lie at a right angle to each other and are connected by at least one fibril. Unlike their closest known relative Tremula longifila, they rotate around their longitudinal axis when swimming and only very rarely glide on surfaces. Screening of a wide range of environmental DNA extractions with lineage-specific PCR primers reveals that Aquavolonida consists of a large radiation of protists, which are most diversified in freshwater planktonic habitats and as yet undetected in marine environments. Earlier-branching lineages in Aquavolonida include less frequently detected organisms from soils and freshwater sediments. The 18S rRNA gene phylogeny suggests that Aquavolonida forms a common evolutionary lineage with tremulids and uncharacterized 'Novel Clade 12', which likely represents one of the deepest lineages in the Rhizaria, separate from Cercozoa (Filosa), Endomyxa, and Retaria.